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ABSTRACT 

f] Carinae (77 Car) is a stellar binary system with a period of 5.54 years. It 
harbors one of the brightest and most massive stars in our galaxy. This paper 
presents spectroscopic evidence for a fast (up to 2,000 km s _1 ) X-ray outflow of 
ionized gas launched from rj Car just before what is believed to be the binary 
periastron (point of smallest binary separation). The appearance of this high- 
velocity component, just as the irregular flares in the X-ray light curve, can not 
be explained by the simple continuous binary wind interaction, adding to the 
intrigue of the rj Car system. 
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- stars: individual: rj Carinae 
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1. INTRODUCTION 



Despite being studied for many yea rs and in all waveba nds, the 77 Carinae system (77 Car), 
best known for its spect acular nebula (IMorse et al. Hl998h that can be tracked back to the 
Great Eruption of 1840 ( Davidson fc Humphreys! Il997l ) is still much of a mystery. It is now 
believed to be a stellar binary system with a period of 5.54 years (jDaminehl Il996l ). The 
primary star is one of the brightest and most massive (~ 120 M@) stars in our galaxy. 
The system is currently in a phase of immense mass ejection. The X-rays detected from 
77 Car can be explained by the collision of the two stellar winds (j Corcoran et al.l l2001al : 
Pittard fc Corcoran! 120021) . not so the 70 day X-ray minimum observed periodically around 
periastron (jCorcoranl 120051 ) . 



Currently, the mass loss rate of the primary excee ds 10 4 M^ vr 1 (jSmith et al.l 120031 ) 



and that of the secondary is an order of magnitude lower (IPittard fc Corcoran 



unusu ally high. The wind of the secondary is much faster (few 1,000 km s 



20o3, but still 



Pittard fc Corcoran 



20021 ) than that of the primary (~ 600 km s 1 and depending on latitude. ISmith et al.ll2003l ). 
The X-ray flux between 2 and 10 keV during the past 12 years is normally at a level of 
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5x10 11 erg s 1 cm 2 , as expected from the co llision of the two stellar winds (IUsovlll992 



Corcoran et al.ll2001al ; iPittard &: Corcoranll2002l ). The X-ray light curve is roughly constant 
throughout most of the 5.54 year orbital period, but rises g radually by a factor o f 3-4 upon 
approach to periastron at which time short flares appear (ICorcoran et al.l 119971 ). before it 
drops sharply and stays in an X-ray low state (< 20% of no rmal brigh t ness) for approx- 
imately 70 days. The full X-ray light curve can be found in ICorcoran! (120051 ) . Although 
absorption by the dense primary wind may play a certain role in attenuating the X-ray 
emission from r\ Car, both spectral and temporal arguments have been put forward to reject 



absorption or an eclipse as the prime reason for this persisting low state (jHamaguchi et al. 



20041 ; ISokerl 120051 ; lAkashi et al.ll2006l ). which is unique to rj Car. Recently, it has been pro- 
posed that if the binary separation during periastron is small enough, the mass ive primary 



wind could smothe r the secondary wind and thus shut down the X-ray emission (jSokerll2005 



Akashi et al.ll2006l ). However, so far there has been no direct evidence for accretion. 



In this work, we wish to add to the temporal aspects of the r\ Car tale a high-resolution 
X-ray spectroscopic dimension. This is achieved by using five deep exposures of r\ Car with 
the High Energy Grating Spectrometer (HETGS) on board the Chandra X- ray Observatory 
(CXO), only the first of which has been published ( Corcoran et al. 2001al lbh . Theoretical 
variations of spectr al line profiles in binary colliding wind systems have been modeled by 



Henley et al.l (120031 ) . The varying line profiles presented here for 77 Car, however, do not 
follow these models, at least not in a straightforward way, as we explain in §3. 



2. Observations & Results 



The log of the Chandra observations is given in Table 1. Two observations (3745 and 
3748) coincide with the short intense flares just before periastron. The exact time of pe- 
riastron is not well determined and depends on the uncertain orbital parameters of the 
system. Nonetheless, for the sake of presentation, we assume throughout this paper a period 
of 5.54 years (2024 days) and set phase zero to be at 2003-06-29. The X-ray low state, thus, 
occurs during phase - 0.035. The data were retrieved from the CXO archive and pro- 
cessed using standard software (CIAO version 3.2.2.). Chandra's HEG and MEG (high and 
medium energy gratings) spectrometers operate simultaneously in the 3.5 - 7.5 A band of 
interest here and their spectra are consistent with each other. Aiming at the highest possible 
kinematic resolution, however, the higher spectral resolving power of the HEG in this band 
(factor of ~ 2) is favored over the slightly higher effective area of the MEG (factor of ~ 2), 
especially since the co-added plus and minus 1st order HEG data are already of sufficiently 
high signal-to- noise ratio (S/N). 
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To reveal the dynamics of the X-ray gas, we transferred each spectral line profile to 
line-of-sight velocity space (v\\ in km s^ 1 ) using the simple non-relativistic Doppler shift 
vu = c(A — A )/A , where A and A represent the observed and rest-frame wavelengths of 
the line, respectively, and c is the speed of light. Negative and positive velocities represent 
approaching and receding gas, respectively. Fig. 1 shows the broad and variable X-ray line 
profiles observed in r] Car as demonstrated on the Si +13 Lya unresolved doublet at 6.18 A. 
As expected, the MEG profile represents a smoothed version of the higher-resolution HEG 
profile. The observed line is much broader than the instrumental line spread function (also 
shown). Fig. 1 also demonstrates how the line profile varies dramatically from assumed 
phase 0= -0.470 (Obs. 632) to 0= -0.028 (Obs. 3745). The early-phase (0= -0.470) line 
is centered around zero velocity with no significant emission beyond 700 km s" 1 (FWHM = 
1000 ± 100 km s _1 ). The line closer to periastron (0= -0.028) features high velocity gas up 
to ~2000 km s _1 . A gaussian fit yields a centroid shift of -550 ± 90 km s _1 , although the 
profile is clearly asymmetric. 

In each observation, the kinematic profiles of all lines are fairly similar. This is demon- 
strated in Fig. 2, where we present four bright spectral lines from Obs. 3745 (0= -0.028). 
The profiles are unambiguously broad and asymmetric with approaching velocities reach- 
ing 2,000 km s _1 , but with no significant receding emission beyond ~ 700 km s^ 1 . In 
order to understand how the line profiles vary between observations and as a function of 
phase we co-added the profiles of nine bright lines in each observation, namely the Lya 
(ls-2p), Hea resonance (ls 2 -ls2p) and Hea forbidden (ls 2 -ls2s) lines of the Si, S, and Ar 
ions. Since the profiles of all these lines in a given observation are roughly similar (e.g., 
Fig. 2), co- adding gives an average profile with improved S/N. Individual Fe lines below 
2 A could not be included in this analysis, despite their high intensity, since the HETGS 
resolving power decreases strongly at these low wavelengths to the extent that the He-like 
Fe line complex is unresolved. The blue wing of the Hea Fe resonance line is consistent 
with, though not resolved as well as, the other line profiles presented here. The resulting 
mean profiles have been normalized (except for Obs. 3747) to facilitate the comparison and 
are presented in Fig. 3. The observed, mean, unnormalized peak intensities for the assumed 
phases: -0.470, -0.130, -0.028, -0.006, and +0.044 are, respectively, 2.0, 3.9, 3.9, 1.9, and 
0.3, x 10 -3 ph s _1 cm~ 2 A -1 . However, any interpretation of the absolute line fluxes in 
terms of the orbital parameters would be risky, as the observed variability is a result of both 
abrupt flaring and varying absorption. 

A clear trend can be seen in Fig. 3. Far from periastron, the line profile is relatively 
symmetrical and narrow. Closer to periastron, bright components of gas moving towards us 
at velocities as high as -2,000 km s _1 start to develop. As the system further approaches 
periastron, the outflow dominates the line profile to the point where the bulk of the emission 
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is clearly blue-shifted. After periastron absorption considerably attenuates the emission lines 
to the low level seen in Fig. 3. Naive gaussian fits to the first two profiles yield centroids at 
-25 ± 25 km s" 1 (for <\>= -0.470) and at -50 ± 20 km s' 1 0= -0.130) both with FWHM = 
700 ± 50 km s _1 . The late-phase asymmetric profiles are non-gaussian. For the mere purpose 
of quantifying the centroid shift, we performed two-gaussian fits, which yield peaks at -120 
± 50 km s" 1 and -800 ± 100 km s" 1 for 0= -0.028 and at -170 ± 50 km s" 1 and -1000 
± 100 km s _1 for <p= -0.006. It is hard to rule out high- velo city receding gas before the X-ray 
minimum, since in agreement with lHamaguchi et al.l (120041 ). we find that the column density 
toward the X-ray source rises from 5xl0 22 cm~ 2 long before periastron to 3xl0 23 cm -2 just 
after it. Emission from the far side of the system, thus, could be strongly absorbed. 



3. Interpretation & Discussion 

The colliding wind binary model, which is generally accepted for explaining the X-ray 
emission of rj Car, has several clear predictions regarding the spectral line profiles expected 
at different orbital phases. The Doppler shifts are the result of shocked gas flowing along the 
contact discontinuity (CD) surface and away from the stagna tion point (SP). If the shock 



opening angle is between 45° —90° ( IPittard Corcoranll2002l ). the broadest line profiles are 



expected when the system is observed perpendicular to the line connecting the two stars, as 
the hottest gas formed around the SP flows directly towards and away from the observer. In 
this case, the lines are centered at zero velocity and are intrinsically symmetrical, but could 
be skewed by absorption of the far (red) emission. On the other hand, when the system 
is observed from behind one of the stars, shifted centroids are expected. If the system is 
observed from behind the weaker (stronger) wind, downstream shocked gas is flowing along 
the CD surface in the general direction of (opposite) the observer, producing blue- (red- 
) shifted centroids. The maximum velocities in this case would be less than those in the 
first (orthogonal) orientation. For a complete and detailed calculation of line profiles from 
colliding winds see Henley et al. (2003). In the following, however, we argue that the current 
line profiles observed at the different phases of the rj Car orbit are in contradiction, even 
qualitatively, with the simple colliding wind scenario. 



If the major axis of the rj Car binary is perpendicular to our line of sight (ISmith et al. 



20041 ). the broadest profiles are expected at apastron (see above), but that is when the 
observed profiles are narrowest (0 = —0.470). Also, with this geometry, the secondary 
would have t o pas s in front of the primary before periastron (in contrary to the picture of 



Smith et al.l (120041 )) to explain the blue-shifted centroids at <\> = -0.028, -0.006, but that 



would produce symmetrical (blue-shifted) profiles as absorption through the secondary wind 
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is weak, while strongly asymmetrical profiles are observed. In any case, the lower velocities 
observed near apastron appear to rule out this orientation if the colliding wind geometry is 
to produce the observed X-ray line profiles. 

If alternatively, the projection of the line of sight on the orbital plane is more or 
l ess along the major ax is and the viewing angle at apastron is from behind the secondary 
( jCorcoran et al.l l2001al ) . the observed blue-shifts at <fi = —0.028, —0.006 could be ascribed 
to gas flowing from the SP towards the observer and the unobserved red wing of the line to 
absorption by the primary wind. However, according to this scenario, at apastron one would 
expect blue-shifted centroids of at least a few 100 km s -1 due to downstream gas (see above), 
which is not observed (0 = —0. 470). For an observer above the binary plane (i < 90°), as 
seems to be the case in rj Car (jCorcoran et al.l l2001al ). the problem becomes more severe 
(higher blue-shift is expected, but not observed). 

In short, the observed behavior of the X-ray line profiles of r\ Car appear qualitatively 
inconsistent with a naive wind-collision scenario, since the interpretation of the high blue - 
shifts from a post-shock stream can not be made consistent with all phases. iHenleyl (120051 ). 
comparing with line profiles calculated from hydrodynamical simulations, also reached the 
conclusion that the observed profiles are not as expected from simple models of the geometry 
of the wind-wind collision. He suggested that including the Coriolis effect may bring the 
models to better agreement with observations. iKashi fc Sokerl (120071 ) suggest the effect of 
orbital motion, which can account for small centroid shifts, but can not explain the currently 
observed X-ray profiles. Here, we wish to propose an alternative explanation. In particular, 
we note that the asymmetric line profiles are observed during strong peaks in the X-ray 
l ight curve. These peaks too can not be acco unted for by the simple colliding wind model 
(IPittard fc Corcoran! |2002| ; iHenley et al.l 120031 ) . or by the Coriolis force. 



We therefore prefer to interpret the observed high-velocity gas up to -2000 km s" 1 in 
terms of a transient, collimated fast wind, or a jet, ejected from the immediate vicinity of 
the binary system through the interaction of the two stars. The peaks in the light curve hint 
that the outflow may be in the form of blobs. The unresolved Chandra images and the rapid 
variations both indicate that the outflow is restricted to within approximately 2xl0 16 cm 
from the center (assuming 0.5" at rj Car's distance of 2.3 kpc). The appearance of the same 
charge states in the spectra attained during all phases implies that the temperatures of the 
X-ray gas and line-emitting outflow remain in the range of kT = 2-5 keV. Consequently, 
we speculate that the fast component of the outflow consists of gas shocked by the collision 
of the winds (as observed throughout the orbit) that is launched and collimated near peri- 
astron. The widths of the major peaks seen in Figs. 2 and 3 suggest that the outflow is only 
moderately collimated to within ~ 30°. The condensation of shocked gas provides a natural 
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explanation for the rise in X-ray intensity during the short intermittent flares observed in 
the X-ray light curve. A collimated outflow from the secondar y has also been sugg ested 
to explain the enhanced He II A4686 emission before periastron (ISoker fc Beharll2006l ). In- 
terestingly, both proper motion and Doppler shifts h ave been measured for high-velocity 
visibl e-light knots much further out from the center (IWalborn et al.l 1 19781 ; iMeaburn et al. 
19931 ). The op tical knots were eject ed more than a hundred years ago along the minor axis 
of the nebula (IMeaburn et al.lll993l ). We would expect the present outflow to be launched 
along the major axis of the nebula although we have no pertinent information on its present 
direction with respect to the system's geometry. If indeed t hat is the cas e and if the inclina- 
tion angle of the binary plane to our line of sight is i = 42° (jSmithll2002l ). the actual outflow 
velocity would be as high as ~2,700 km s _1 . 

As the massive r\ Car primary is a short-lived star shedding considerable mass and as it 
is now understood that core collapse supernovae can not expel more than a few solar masses, 
implying that higher mass stars must shed most of their mass prior to the explosion, the winds 
and fast outflow of rj Car may be a supernova in the making. At the very least, it highlights 
the next periastron passage on 2009, January 12 as a faithfully scheduled experiment for the 
astrophysically common, but poorly understood phenomenon of jet launching. 
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Table I. CXO/HETG Observation Log 



CXO archive 


Start Time 


Exposure 


Assumed Orbital 


2-10 keV Flux 


ID a 




(ks) 


Phase 0) b 


(10 -10 erg s -1 cm -2 ) 


632 


2000-11-19 02:46:40 


90.69 


-0.470 


0.50 


3749 


2002-10-16 08:08:49 


93.96 


-0.130 


0.98 


3745 


2003-05-02 11:56:16 


97.29 


-0.028 


2.2 


3748 


2003-06-16 05:35:28 


100.1 


-0.006 


0.97 


3747 


2003-09-26 22:45:53 


72.16 


+0.044 


0.48 



a One additional observation was carried out during the X-ray minimum (ID 3746 at 2003- 
07-20 01:46:22) for which no useful spectrum could be obtained. 

b Phase is approximated based on the assumption of a 5.54 year (2024 day) orbit, where 
phase zero is set at 2003-06-29 and the X-ray low state occurs during phase - 0.035. 
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Line-of-sight velocity (km s" 1 ) 



Fig. 1. — Velocity profiles of the Si +13 Lya unresolved doublet (6.180 and 6.186 A). Upper 
panel shows consistent HEG (higher spectral resolution) and MEG (smoother) data from 
Obs. 3745 (<f)= -0.028). The profile is clearly broadened up to 2000 km s _1 , much beyond 
the instrument line spread functions (dashed lines). Bottom panel shows the HEG profile of 
Obs. 632 (0= -0.470, dashed line) and Obs. 3745 (0= -0.028). The scaled up line from 
the early phase is rather symmetrical and much narrower, showing no line emission beyond 
700 km s~ l . 
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Fig. 2. — Profiles of four bright spectral lines from Obs. 3745 (assumed phase -0.028). The 
consistent asymmetric profiles with blueshifts up to ~ —2, 000 km s _1 are clearly seen. 
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Fig. 3. — Mean, normalized, velocity profiles constructed from nine different spectral lines 
for the five observations of i] Car. Only the profile of Obs. 3747 (assumed phase +0.044) 
during which rj Car was strongly absorbed is not normalized and plotted at its correct scale 
with respect to Obs. 632 (assumed phase -0.470). Typical la errors on each data point are 
10 - 15 % (top right hand side). 



